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Abstract 

By investigating model-independent bounds for exotic options in financial mathe- 
matics, a martingale version of the Monge-Kantorovich mass transport problem was 
introduced in [3, 24]. In this paper, we extend the one-dimensional Brenier's theorem 
to the present martingale version. We provide the explicit martingale optimal transfer- 
ence plans for a remarkable class of coupling functions corresponding to the lower and 
upper bounds. These explicit extremal probability measures coincide with the unique 
left and right monotone martingale transference plans, which were introduced in [4] by 
suitable adaptation of the notion of cyclic monotonicity. Instead, our approach relies 
heavily on the (weak) duality result stated in [3], and provides, as a by-product, an 
explicit expression for the corresponding optimal semi-static hedging strategies. We 
finally provide an extension to the multiple marginals case. 

1 Introduction 

Since the seminal paper of Hobson [29], an important literature has developed on the topic 
of robust or model-free superhedging of some path dependent derivative security with payoff 
£, given the observation of the stochastic process of some underlying financial asset, together 
with a class of derivatives. See [7, 11, 12, 13, 14, 15, 16, 18, 19, 31, 33, 39] and the survey 
papers of Obloj [40] and Hobson [30] . In continuous-time models, these papers mainly focus 
on derivatives whose payoff £ is stable under time change. Then, the key-observation was 
that, in the idealized context where all T— maturity European calls and puts, with all possible 
strikes, are available for trading, model-free superhedging cost of £ is closely related to the 
Skorohod Embedding problem. Indeed, the market prices of all T— maturity European calls 
and puts with all possible strikes allow to recover the marginal distribution of the underlying 
asset price at time T. 

Recently, this problem has been addressed via a new connection to the theory of optimal 
transportation, see [3, 24, 27, 1, 2, 20, 21]. Our interest in this paper is on the formulation of 
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a Brenier Theorem in the present martingale context. We recall that the Brenier Theorem 
in the standard optimal transportation theory states that the optimal coupling measure is 
the gradient of some convex function which identifies in the one-dimensional case to the so- 
called Frechet-Hoeffding coupling [6]. A remarkable feature is that this coupling is optimal 
for the class of coupling cost functions satisfying the so-called Spence-Mirrlees condition. 

We first consider the one-period model. Denote by X, Y the prices of some underlying 
asset at the future maturities and 1, respectively. Then, the possibility of dynamic trading 
implies that the no-arbitrage condition is equivalent to the non-emptyness of the set A4 2 
of all joint measures P on R + x R + satisfying the martingale condition E P [Y|X] = X. 
The model-free subhedging and superhedging costs of some derivative security with payoff 
c(X,Y), given the marginal distributions X ~ ji and Y ~ z/, is essentially reduced to the 
martingale transportation problems: 

inf E p [c{X,Y)} and sup E F [c{X,Y)], 
PeM 2 0^) PeM 2 (M,^) 

where M2(/i, v) is the collection of all probability measures P G M.2 such that X ~p /z, 
Y ~p v. Our main objective is to characterize the optimal coupling measures which solve 
the above problems. This provides some remarkable extremal points of the convex (and 
weakly compact) set M 2 (/i, v). In the absence of marginal restrictions, Jacod and Yor [35] 
(see also Jacod and Shiryaev [34], Dubins and Schwarz [22], for the discrete-time setting) 
proved that a martingale measure P G M.2 is extremal if and only if P-local martingales 
admit a predictable representation. In the present one-period model, such extremal points 
of M.2 consist of binomial models. For a specific class of coupling functions c, the extremal 
points of the corresponding martingale transportation problem turn out to be of the same 
nature, and our main contribution in this paper is to provide an explicit characterization. 

Our starting point is a paper by Hobson and Neuberger [32] who considered the specific 
case of the coupling function c(x,y) := \x — y\, and provided a complete explicit solution 
of the optimal coupling measure and the corresponding optimal semi-static strategy. In a 
recent paper, Beiglbock and Juillet [4] address the problem from the viewpoint of optimal 
transportation. By a convenient extension of the notion of cyclic comonotonicity, [4] in- 
troduce the notion of left-monotone transference plan. They also introduce the notion of 
left-curtain as a left-monotone transference plan concentrated on the graph of a binomial 
map. The remarkable result of [4] is the existence and uniqueness of the left-monotone 
transference plan which is indeed a left-curtain, together with the optimality of this joint 
probability measure for some specific class Cbj of coupling payoffs c(x,y). Notice that the 
coupling measure of [32] is not a left-curtain, and Cbj does not contain the coupling payoff 
\x - y\. 

As a main first contribution, we provide an explicit description of the left-curtain P* of [4] . 
Then, by using the weak duality inequality, 

- we provide a larger class C D Cbj of payoff functions for which P* is optimal, 

- we identify explicitly the solution of the dual problem which consists of the optimal 
semi-static superhedging strategy, 
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- as a by-product, the strong duality holds true. 

Our class C is the collection of all smooth functions c : R x R — > R, with linear growth, 
such that c xyy > 0. We argue that this is essentially the natural class for our martingale 
version of the Brenier Theorem. 

We next explore the multiple marginals extension of our result. In the context of the 
finite discrete-time model, we provide a direct extension of our result which applies to the 
context of the discrete monitored variance swap. This answers the open question of optimal 
model-free upper and lower bounds for this derivative security. 

The paper is organized as follows. Section 2 provides a quick review of the Brenier Theorem 
in the standard one-dimensional optimal transportation problem. The martingale version of 
the Brenier Theorem is reported in Section 3. We next report our extensions to the multiple 
marginals case in Section 4. Finally, Section 5 contains the proofs of our main results. 



2 The Brenier Theorem in One-dimensional Optimal 
Transportation 

2.1 The two- marginals optimal transportation problem 

Let X, Y be two scalar random variables denoting the prices of two financial assets at 
some future maturity T. The pair (X, Y) takes values in IR 2 , and its distribution is defined 
by some probability measure P G Pr2, the set of all probability measures on IR 2 . For the 
purpose of the present financial application, the measures have support on R+. For the sake 
of generality, we consider however the general case. 

We assume that T— maturity European call options, on each asset and with all possible 
strikes, are available for trading at exogenously given market prices. Then, it follows from 
Breeden and Litzenberger [5] that the marginal distributions of X and Y are completely 
determined by the second derivative of the corresponding (convex) call price function with 
respect to the strike. We shall denote by fx and v the implied marginal distributions of X 
and Y, respectively, £ M ,r M , l v \r v the left and right endpoints of their supports, and F^, F v 
the corresponding cumulative distribution functions. 

By definition of the problem the probability measures \i and v have finite first moment: 

J \x\fj,(dx) + J \y\v(dy) < oo, (2.1) 

and although the supports of ji and v could be restricted to the non-negative real line for 
the financial application, we shall consider the more general case where /i and v lie in Vm., 
the collection of all probability measures on IR. 

We consider a derivative security defined by the payoff c(X, Y) at maturity T, for some 
upper semicontinuous function c : IR 2 — > IR with linear growth. This condition could be 
replaced by 

c(x,y) < <p(x) +tjj(y) for some cp,tp:R — >R, cp + e L 1 ^), ^ + e h\iy). (2.2) 
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The model-independent upper bound for this payoff, consistent with vanilla option prices of 
maturity T, can then be framed as a Monge-Kantorovich (in short MK) optimal transport 
problem: 

po^, v ) := sup E p [c(X, Y)] where P 2 (/i, u) := {P G V R 2 : X ~ P fi and Y ~ P z/}, 

where, for the sake of simplicity, we have assumed a zero interest rate. This can easily be 
relaxed by considering the forwards of X and Y. Notice that c(X, Y) is measurable by the 
upper semicontinuity condition on c, and is integrable by the linear growth condition on c 
together with (2.1) (or the condition (2.2)). 

In the original optimal transportation problem as formulated by Monge, the above maxi- 
mization problem was restricted to the following subclass of measures. 

Definition 2.1. A probability measure P G ^(/-f, v) is called a transference map ifP(dx, dy) : = 
fi(dx)5{T(x)}(dy) , for some measurable map T : K. — > R. 

The dual problem associated to the MK optimal transportation problem is defined by : 

D° 2 (n,v) := inf Ufa) + vty)}, 

where, denoting ip © tp(x, y) := <f(x) + ip(y): 

V\ := {((p,ip) : e V-(n),il> + G V{v) and © ^ > c}. 

and with /x(</?) := / </?d/x, z/(-^) := J -^rfz/. 

The dual problem D^ifi, v) is the cheapest superhedging strategy of the derivative security 
c(X, Y) using the market instruments consisting of T— maturity European calls and puts 
with all possible strikes. The weak duality inequality 

P 2 <W) < D° 2 (fi,u) 

is immediate. For an upper semicontinuous payoff function c, equality holds and an optimal 
probability measure P* for the MK problem P 2 ° exists, see e.g. Villani [43]. 

Our main interest of this paper is the following one-dimensional version of a result estab- 
lished by Brenier [6], which provides an interesting characterization of P* in terms of the 
so-called Frechet-Hoeffding pushing forward /j, to v, defined by the map 

T* := F- 1 o F„ (2.3) 

where F~ l is the right-continuous inverse of F v : 

F~\t) := inf{y : F v (y) > x}. 

In particular, the following result relates the MK optimal transportation problem P® to the 
original Monge mass transportation problem for a remarkable class of couplings c. This result 
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is more general, in particular the set of measures Ft induced by a map T pushing forward \i 
to v is dense in V^p. whenever \x is atomless and we consider compact subsects of R 2 . For the 
purpose of our financial interpretation, this result characterizes the structure of the worst 
case financial market that the derivative security hedger may face, and characterizes the 
optimal hedging strategies by the functions </?* and tp* defined up to an irrelevant constant 
by 

<p*(x):=c(x,T*(x))-rl>.oT*(x), ^(y) := c y (T-\y),y) , x,y ER. (2.4) 

Theorem 2.2. (see e.g. [43], Theorem 2.44) Let c be upper semicontinuous with linear 
growth. Assume that the partial derivative c xy exists and satisfies the Spence-Mirrlees con- 
dition c xy > 0. Assume further that /i has no atoms, G and ip+ G L 1 (z/). Then 

(i) P°(ti,u) = D° 2 (fi,is) = fc(x,T m (x))ii(dx), 

(ii) (</?*, ■0*) £ 12 2, and is a solution of the dual problem D%, 

(hi) F*(dx,dy) := /i(dx)5T*( x )(dy) is a solution of the MK optimal transportation problem 
P®, and is the unique optimal transference map. 

Proof. We provide the proof for completeness, as our main result in this paper will be 
an adaptation of the subsequent argument. First, it is clear that P* G V(fi,u). Then 
E p *[c(X, Y)} < P 2 °(/i, v). We now prove that 

(y?*,^)GX?° and //(<£>*) + i^*) = E p * [c(X, Y)]. (2.5) 

In view of the weak duality -P 2 °(/i, u ) ^ -^(/^ u), this would imply that -P 2 °(/i, v ) = -^(/^ v ) 
and that P* and (</?*, ip*) are solutions of P 2 (A 4 ? z/ ) an d ^(fi, v), respectively. 
Under our assumption that </?* G L 1 (/i), ^+ G L 1 (z/), notice that (2.5) is equivalent to: 

= H°(x,T*(x)) = minH°(x,y), where H° :=</?*© V* - c. 

3/GK 

The first-order condition for the last minimization problem provides the expression of ip^ in 
(2.4), and the expression of ip* follows from the first equality. Since 

rT«Hy) 

Hy{x,y) = c y {T; l (y),y) - c y (x,y) = / c xy {£,y)d£, 

J X 

it follows from the Spence-Mirrlees condition that T* (x) is the unique solution of the first- 
order condition. Finally, we compute that Hy y {x,T*(xj) = c xy (x, T*(xj) /Tl(x) > by the 
Spence-Mirrlees condition, where the derivatives are in the sense of distributions. Hence 
T*(x) is the unique global minimizer of H(x, .). O 

We observe that we may also formulate sufficient conditions on the coupling c so as to 
guarantee that the integrability conditions y9+ G L 1 (//),'0+ G L 1 (z/) hold true. See [43], 
Theorem 2.44. 
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Remark 2.3 (Mirror coupling: anti-monotone rearrangement map), (i) Suppose that the 
coupling function c satisfies c xy < 0. Then, the upper bound P^fi,^) is attained by the 
anti-monotone rearrangement map 

W*(dx,dy) := n(dx)S {Tt{x)} (dy), where T m (x) := F' 1 o (l - F^(-x)). 

To see this, it suffices to rewrite the optimal transportation problem equivalently with modified 
inputs: 

c(x,y) := c(-x,y), Jl(x) := n((-x, oo)), V := v, 
so that c satisfies the Spence-Mirrlees condition c xy > 0. 

(ii) Under the Spence-Mirrlees condition c xy > 0, the lower bound problem is explicitly solved 
by the anti-monotone rearrangement. Indeed, it follows from the first part (i) of the present 
remark that: 

inf E p [c(X,Y)]=- sup E ¥ [-c(X,Y)] = -E ¥ *[- c(X,Y)] = [c(x,T*(x))fi(dx). 

PeTM^.") "' Pe7> 2 (M,") " J 

Remark 2.4. The Spence-Mirrlees condition is a natural requirement in the optimal trans- 
portation setting in the following sense. The optimization problem is not affected by the 
modification of the coupling function from c to c := c + a © b for any a G and 
b G L 1 (z/). Since c xy = c xy , it follows that the Spence-Mirrlees condition is stable for the 
above transformation of the coupling function. 

Example 2.5 (Basket option). Let c(x,y) = (x + y — k) + , for some k G K. (see [17, 38] 
for multi-asset basket options). The result of Theorem 2.2 applies to this example as well, 
as it is shown in [43] Chapter 2 that the regularity condition c G C 1,1 is not needed. The 
upper bound is attained by the Frechet-Hoeffding transference map T* := F~ l o F^, and the 
optimal hedging strategy is: 

^*(y) = (y - y) + , v*( x ) = ( T *(X> +x-k) + - (t*(x) - y) + , 

where y is defined by T*(k — y) — y. 

2.2 The multi-marginals optimal transportation problem 

The previous results have been extended to the n— marginals optimal transportation problem 
by Gangbo and Swi§ch [25], Carlier [9], and Pass [41]. Let X = (X 1; . . . ,X n ) be a random 
variable with values in W 1 , representing the prices at some fixed time horizon of n financial 
assets, and consider some upper semicontinuous payoff function c : W 1 — > R with linear 
growth. 

Let /ii, ...,/t n e Pr be the corresponding marginal distributions, and \i := (/xi, . . . , /i n ). 
The upper bound market price on the derivative security with a payoff function c is defined 
by the optimal transportation problem: 

P n °(/i) := sup E p [c(X)] , where V n (fi) := {P G TV : X t ~ P fa, 1 < i < n}. (2.6) 
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Then, under convenient conditions on the coupling function c (see Pass [41] for the most 
general ones), there exists a solution P* to the MK optimal transportation problem -P°(/f) 
which is the unique optimal transference map defined by T*, i — 2, . . . , n: 

n 

F*(dx 1 , dx n ) = Hi{dxi) Y[ S Ti{xi) (dxi), where T\ = F' 1 oF W) % = 2, . . . , n. 

i=2 

The optimal upper bound is then given by 

3 The Two- Marginals Martingale Transport Problem: 
Main Results 

The main objective of this paper is to obtain a version of the Brenier theorem for the 
martingale transportation problem introduced by Beiglbock, Henry-Labordere and Penkner 
[3] and Galichon, Henry-Labordere and Touzi [24]. A result in this direction was first 
obtained by Hobson and Neuberger [32] and by Beiglbock and Juillet [4]. In contrast with the 
last reference, our result is an explicit extension of the Frechet-Hoeffding optimal coupling. 
We outline in Sections 3.6 and 3.7 the main differences with [4, 32]. 

3.1 Problem formulation 

In the context of the financial motivation of Subsection 2.1, we interpret the pair of random 
variables X, Y as the prices of the same financial asset at dates ti and t 2) respectively, with 
t\ <t 2 . Then, the no-arbitrage condition states that the price process of the tradable asset 
is a martingale under the pricing and hedging probability measure. We therefore restrict 
the set of probability measures to: 

M 2 (n,v) := {F eV 2 (fi,v) :E P [Y\X]= X}. 

where /i, v have finite first moment as in (2.1). This set of probability measures is clearly 
convex, and the martingale condition implies that t < < r M < r v ' . Throughout this 
paper, we shall denote 

SF := F v - Fp. 

By a classical result of Strassen [42], M. 2 {^-, v) is non-empty if and only if \x < v in sense of 
convex ordering, i.e. 

(i) /i, v have the same mean: J ^d5F(^) = 0, 
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(ii) and /(£ - k) + /j(d£) < /(£ - k) + v(d£), for all fc G R. This condition can also be 
expressed as: 

A 5F{£)d£<0 or, equivalently, j 5F{£)d£>0, for all k G R, (3.1) 

«/[fc,oo) J [— oo,fc) 

where the last equivalence follows from the first property (i). 

For completeness, we provide in Section 6 some examples of probability measures in M. 2 (/i, v) 
which are commonly using by practitioners in quantitative finance. 

Let c : R 2 — > R be an upper semicontinuous function with linear growth (or the condition 
(2.2)), representing the payoff of a derivative security. In the present context, the model- 
independent upper bound for the price of the claim can be formulated as the following 
martingale optimal transportation problem: 

P 2 (fi,v) := sup E ¥ [c(X,Y)], (3.2) 

Remark 3.1. When /i and v have finite second moment, notice that E P [(X — Y) 2 } = 
-E P [X 2 ] + E F [Y 2 ] = J £ 2 d5F(£) for all P G M(ji,v). Then, the quadratic case, which 
is the typical example of coupling in the optimal transportation theory, is irrelevant in the 
present martingale version. 

We finally report the Kantorovich dual in the present martingale transport problem. Be- 
cause of the possibility of dynamic trading the financial asset between times t\ and t 2 , the 
set of dual variables is defined by: 

V 2 := {(</?, ■?/>, h) : (p + E G L 1 (z/), h E L°, and (p ® i> + h® > c) , (3.3) 

where ip © ip(x,y) := (p(x) + if)(y), and h®(x,y) := h(x)(y — x). The dual problem is: 

D 2 (»,v) := inf {fifoO + i/ty)}, (3.4) 

(tp,il>,h)eV2 

and can be interpreted as the cheapest superhedging strategy of the derivative c(X, Y) by 
dynamic trading on the underlying asset, and static trading on the European options with 
maturities t\ and t 2 . Since c has linear growth and //, v have finite first-order moments, the 
weak duality inequality: 

Pi(p,v) < D 2 {n,v) (3.5) 

follows immediately from the definition of both problems. Under suitable conditions on c, [3] 
proved the strong duality result (i.e. equality holds), and showed the existence of a maximizer 
P* G A4 2 (fi, u) for the martingale transportation problem P 2 (/i, v). However, existence does 
not hold in general for the dual problem D 2 (/i, v). An example of non-existence is provided 
in [3]. 
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3.2 Preliminaries 



Our objective in this section is to provide explicitly the left-monotone martingale transport 
plan, as introduced by Beiglbock and Juillet [4]. 

Definition 3.2. We say thatF G A^A*, v) is left-monotone (resp. right-monotone) if there 
exists a Borel set Y C R x R snc/i £/ia£ P[(X, F) G T] = 1, and for all (x, j/i), (x, y 2 ), (x 1 , y') G 
T with x < x 1 (resp. x > x' ), it must hold that y' £ (2/1,2/2)- 

Our main results hold for probability measures /i, v satisfying the following restriction. 

Assumption 3.3. The probability measures /i and v have finite first moments, n < v in 
convex order, and \i has no atoms. 

Under this assumption, Theorem 1.5 and Corollary 1.6 of [4] state that there exists a 
unique left-monotone martingale transport plan P* G M. 2 ^,v\ and that the graph of P* 
is concentrated on two maps T d ,T u : R — > R, T d (x) < x < T u (x) for all x G R, i.e. 
P*[y = T d (X)]+F*[Y = T U (X)} = 1. 

Remark 3.4. The condition that is continuous in Assumption 3.3 implies that 5F is 
upper- semicontinuous , and therefore the local suprema of 5F are attained by maximizers. 

For our construction, we introduce the functions: 

g(x,y) := F~ x (F^x) + 5F(y)) , x,y eR, (3.6) 

where F~ l denotes the right-continuous inverse of F v , with F^ 1 = 00 on (1, 00) and F^ 1 = 
—00 on (—oo,0). We also define for a measurable subset A G £>r such that 5F is increasing 
on A: 

G A (t,x) := [ SdF v (£) - f ZdF^O + / {g(x,0 - C)d5F(0, t<xeR. 

J (-oo,F v 1 oF ll (x)] J -00 JAn(-oo,t] 

(3.7) 

In the last integral, notice that g(x,£,) — C, > 0, so that by the increase of 5F on A, the 
integral has a well-defined value in (—00, 00]. It will be made clear in Section 5.1 that these 
functions appear naturally when one imposes that P* G A^A*, v). 

Notice that G A is right-continuous in t, and G A (— 00, 00) = 0, a consequence of the fact 
that n and v have the same mean. Our construction uses the following preliminary result, 
which needs the additional notation: 

B := {x <E M. : 5F increasing to the right of x}, x := inf B , 

where we say that a function <fi is increasing (resp. decreasing) to the right of x if for all 
E > 0, there exists e G (0, e ) such that <f>(x + e) > <p{x) (resp. <f>(x + e) < 4>{x)). 
Observe that x = 00 if and only if (j, — v. 
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Lemma 3.5. Assume x < oo, let m e R be a local maximizer of 5F, and consider a Borel 
subset A C (xo,m] D B . Denote A m := (xo,m] \ A, and assume that f Am d(p(5F) > for 
any non- decreasing function <fi. Then, there exists a unique scalar t A (x,m) such that, for all 
x > m with SF(x) < 5F{m), 

t A (x,m)eA, and G A (t A (x,m)-, x) < < G A (t A (x,m), x) . 

Moreover, x(m) := inf{x > m : g{x, t A (x, raj) < x} satisfies, whenever x(m) < oo, 

5F{t A (x(m),m)) < 5F{x(m)) < 5F(t A (x(m)-,m)) (3.8) 
and 5F is strictly increasing on a right neighborhood of x(m). 

The proof of this lemma is reported in Subsection 5.1. 
3.3 Explicit construction 

Our explicit construction requires an additional condition on 5F. Let ~M(5F) denote the 
collection of all points m such that SF is nondecreasing to the left of m, and decreasing to 
the right of m: 

M(5F) := {m : SF'(m—) < 0, and 5F decreasing to the right of m}, (3.9) 

where we recall that SF decreasing to the right of m means that any right neighborhood of 
m contains a point w! such that 5F{m') < 5F(m). Our general characterization of the left 
monotone transference plan will be obtained in Theorem 3.11 as a limit of explicit monotone 
transference plans corresponding to an approximating sequence satisfying the following no 
right accumulation requirement. 

Assumption 3.6. M(5F) U {x } has no right accumulation point. 
Lemma 3.7. Under Assumption 3.6, the set M(<5F) is countable. 
Proof. Under Assumption 3.6, we have ~M(5F) = \J ne ?q'M n (5F), where 

M n (<5F) := {m : 5F strictly decreasing on (m,m H — ], and SF'(m—) < 0}. 

Then the required result follows from the fact that ~M n (5F) is countable. O 

We are now ready for the construction of the left-monotone transference map P*. We first 
initialize the construction in Step 0, and continue an iterative construction in the subsequent 
steps. 

Step 0: If xq = —oo, we move to Step 1 of the construction. Otherwise, define: 

T d (x) = T u (x) = x for x < x . (3.10) 
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If xo = oo, i.e. /i — v, then this completes the construction of (T d ,T u ). Otherwise, we 
continue with the following step. 

Step 1: By Assumption 3.6, the function 5F increases at the right of Xq. Consider the 
first point of decrease of SF (see Remark 3.4): 

mi := inf {to > x : 5F decreasing on [m,m + e) for some e > 0}, 

By the stochastic dominance /i ■< v, see (3.1), it follows that SF is nondecreasing on (x , mi], 
5F(mi) > 0, and 5F is strictly increasing on the set 

Ai := (x Q ,mi] n B . 

We have A™ 1 = and we are then in the context of application of Lemma 3.5 with (to, A) = 
(mi,Ai). Denoting xi := x(mi), we define the maps T d ,T u on (— oo,Xi): 

T d (x) = T u (x) = x for x <x<mi, 
T d {x):=t^{x,mi), T u (x) := g(x, T d (x)) for mi < x < x ± . l ' } 

If xi = oo, this completes the construction. See Figure 1 below for such an example. 
Otherwise, it follows from Lemma 3.5 that 5F is strictly increasing at the right of x\, 
whenever xi < oo. In this case, we continue the construction denoting: 

Bi := B \ {T d ([mi,xi)) U [mi,T u (xi))} = B \ (T d (xi), T u (xi)), 

where the last equality follows from the fact that T d is decreasing and T d (x) < x, see Remark 
3.8. 

Step 2: The construction of this step falls in the more general Step % below, and is provided 
here for the convenience of the reader. 

Since /i ■< v in (3.1), it follows that the set of local maximizers after x\ is not empty. 
Recall Assumption 3.6, and let: 

m 2 := inf {m > x\ : 5F decreasing on [to, to + e) for some e > 0}, 
A 2 := (x 1 ,m 2 ]nB 1 = (x , T d (xi)) U (x u m 2 ], 



\rr12 



so that 5F is nondecreasing on [xi,TO2], and strictly increasing on A 2 . Moreover A r 2 
[T d (xi),xi] and, since 5F(T d (xi)) < SF(xi) by (3.8), we see that /^m 2 d(f)(8F) > for all 
nondecreasing function 0. 

Then, we may apply Lemma 3.5 with (to, A) = (m 2 ,A 2 ). Denoting x 2 := x(m 2 ), we may 
define the maps T d ,T u on [xi,x 2 ): 

T d (x) = T u (x) = x for xi < x < m 2 , , . 

T d (x) := t M (x, to 2 ), T u (x) := g(x,T d (x)) for m 2 < x < x 2 . 

If x 2 = oo, this completes the construction. See Figure 2 below for such an example. 
Otherwise, it follows from Lemma 3.5 that SF is strictly increasing at the right of x 2 , 
whenever x 2 < oo. In this case, we continue the construction denoting: 

B 2 :=B 1 \{T d ([m 2 ,x 2 )) U [to 2 , T u (x 2 )) }. 
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Step i: Suppose that (T d , T u ) are defined on (—00, xi) for some Xi with SF strictly increas- 
ing at the right of Xi, and let a subset Bi := B^i \ {T d ([rrii, x^) U [m*, Xj]} C -Bo be given so 
that by definition, we have 

G Ai (T d (xi),Xi) > 0. (3.13) 

and Ai is obtained iteratively from the previous steps as: 

A,i = (x ,mi\ \ [Uj<i {^([mj,^)) U [m^, T^)) }] . 

Since fi ^ u in (3.1), it follows that the set of local maximizers after Xi is not empty. Recall 
Assumption 3.6, and let: 

m i+ i := inf jm > X; L : SF decreasing on [m, m + e) for some e > 0}, 

and 

A+i : = (xo.mi+ijnB, = (x ,m i+ i] \ [Uj<i {T d ([m j: Xj)) U [m^, T^-)) }] , 

so that SF is strictly increasing on Ai+\. 

We observe that T d {xi) [T d (xj),T u (xj)\ for any j < i, which expresses that our con- 
struction provides the left-monotone martingale transport plan, see Definition 3.2. Since 
8F(T d (xi)) < SF(xi) by (3.8), we have also that J A m i+i d<p{5F) > for all nondecreasing 

i+l 

function <fi. 

We have thus verified that the conditions of Lemma 3.5 are satisfied by the pair {m i+1 , A i+1 ), 
and we may then define the maps T d ,T u on [xi,x i+ i) by: 

T d (x) = T u (x) =x for Xi < x < m i+1 , 
T d (x) := t Ai+1 (x,m i+1 ), T u (x) := g(x,T d (x)) for m i+1 < x < x i+1 := x(m i+1 ). 

If Xi + \ = 00, the construction is complete. Otherwise, it follows from Lemma 3.5 that SF is 
strictly increasing at the right of Xj+i, whenever Xi + \ < 00. In this case, we also update: 

B l+l := Bi \ {T d ([m i+ i,a; i+ i)) U [m i+1 , T u (x i+1 ))}, 

and we continue with an additional step. 

Case of accumulation : It may happen that the increasing sequence (m^ converges 
to some m l < 00. Then, as the number of steps i tends to infinity, the above construction 
defines the maps (T d ,T u ) on (— 00, m 1 ). 

In this case, under Assumption 3.6 which excludes any right accumulation of local maxima, 
we may start again the construction exactly as in Step i, with m i+1 = m\. After possibly i 
steps, this defines (mj,xj)j<j which either meets the requirement x\ = 00, or accumulates. 
Recall that the set M(5) of (3.9) is countable under our Assumption 3.6. Since the set of 
possible accumulation points m k is a subset of ~M(5F), it is at most countable. Then, by 
transfinite induction, 

we relabel the sequence {m k ^x^)j^ as a new sequence that we rename (m*, Xj)j>o. 
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Remark 3.8 (Some properties of Td). From the above construction ofTd, we see that 

(i) Td is right- continuous. Moreover, on each interval (m,i,Xi), it is non-increasing and flat 
if and only if it reaches an atom of F v . 

(ii) In general, the restriction ofTd to Uj>o(mj, Xj) fails to be non- decreasing. However, for 
i 7^ j, we have Td((rrii,Xi)) DTd((mj,Xj)) = 0. Consequently, the right- continuous inverse 
T^ 1 ofTd is well defined. 

(iii) Let I = (a,b) C T d ([mi,Xi\) be such that 5F is flat on I, and 5F increases at the right 
of b and at the left of a. Then, whenever T d reaches the right endpoint b, it jumps from b to 
a, i.e. AT d (T^(b)) = a - b. 

(iv) Let x be such that Td(x) = T u (x) = x. Then, {x' ^ x : T u {x') = x} = 0, {x' ^ x : 
T d (x') = x} 7^ 0, and reduces to a single point set if AF„(x) = 0. Otherwise, if x is an atom 
of F v , the last set has a positive measure under F M . 

Remark 3.9 (Some properties of T u ). From the above construction ofT u , we see that 

(i) T u ([m h Xi\) C [m h Xi\, and T u (x) > x for x e {m h Xi) for alii. 

(ii) T u is right continuous with discontinuity points {x : AF u (T d (x)) > 0}, recall that F^ is 
continuous. 

(iii) T u is nondecreasing, and strictly increasing on the support of fi. The last property will 
be clear from Theorem 3.10 (ii) below, and implies that the right- continuous inverse T^ 1 of 
T u is well-defined. 

3.4 The left-monotone martingale transport plan 

The last construction provides our martingale version of the Frechet-Hoeffding coupling: 

%(x, dy) := l D (x)5 {x} (dy) + l D c(x) [q(x)5 {Tu(x)} (dy) + (1 - q)(x)5 {Td{x)} (dy)] , (3.15) 
where x_i = — oo, m Q := x , 

D := U i > (a; i _i,m i ] and q(x) := - X . (3.16) 

Observe that T d (x) < x < T u (x) from our previous construction. Therefore, q takes values 
in [0,1]. 

Theorem 3.10. Let Assumptions 3.3 and 3.6 hold true. Then, 

(i) the probability measure P*(ofo;, dy) := n{dx)T*{x, dy) is the unique left-monotone transport 
plan in M. 2 {n, v); 

(ii) moreover T u and T d solve the following ODEs: 

d(5F o Td) = —(1 — q)dF^, d(F u o T u ) = qdF^ whenever x e [m,i,Xi) and T d (x) e int(Aj). 

The proof is reported in Section 5.1. The next result characterizes the left-monotone trans- 
ference map in the case where SF does not satisfy Assumption 3.6. 
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Theorem 3.11. Let Assumption 3.3 hold true, and let (/i n , v) n >\ C Vr be such that \i n — > fi 
and v n — > v , weakly, and (/i n , u n ) satisfies Assumptions 3.3 and 3.6. For all n > 1, define 
the corresponding T™ as in (3.15), and the corresponding F™(dx,dy) := fi n (dx)T™(x,dy) . 
Then P™ converges weakly towards the unique left-monotone transference map. 

Proof. By following the proof of Proposition 2.4 of [3], it follows from Lemma 4.4 p56 in 
Villani [43] that the sequence (P") n >i is weakly compact. Then, after possibly passing to 
subsequence, P" — > P, weakly, for some P G -M(/i, v). To prove the required result, we 
shall prove that P is a left-monotone transference map; then, from the uniqueness result of 
Theorem 1.5 in [4], we may deduce that P does not depend on the chosen subsequence. 

Assume to the contrary that P is not left-monotone. Then there exists a support T of P 
such that 

(x,y d ), (x,y u ), (x',y) G f , y d < y u , x > x, and y G (y d ,y u )- (3.17) 

To obtain the required contradiction, we prove below that there exist sequences (x n ,y d l ) n , 
(x n ,y™) n , (x' n ,y' n ) n in a support of P™ such that (x n ,x' n ) — > (x,x'), and (y d ,yu,y' n ) — y 
(y d , y u , y')- By the left-monotonicity of P™ for all n, we have y' n ^ (y^, y™), and we obtain by 
sending n — > oo that y' £ (yd,Uu), contradicting (3.17). 

We finally prove that if (x,y) G f, then there exists a sequence (x n ,y n ) and a support of 
P n such that (x n ,y n ) — > (x,y). For an arbitrary e > 0, let ip a continuous function with 
support in B e (x,y), the open ball centered at (x,y) with radius e. Then, it follows from 
the weak convergence of P™ towards P that E p * [<p(X, Y)] — > EF[(p(X, Y)], and the required 
result follows from the arbitrariness of e > 0. O 

We conclude this subsection by the following remarkable property of T d . 

Proposition 3.12. Let Assumptions 3.3 and 3.6 hold true. Let i > 1 be such that 5F is 
not flat at the left ofrrii. Then T d (rrii-\-) = mi. If in addition F^jF^ are twice differentiable 
near rrii, then: 

T' d (mi+) = -1/2 and T' d \m,A = +oo. 

Proof By construction, we have T d (rrii+) = rrii. Denoting e := x—T d (x), f^ := F' , f v := F' v , 
8f '■= fv — ftii and recalling that g[x,x) = x, we see by direct calculation that 

SfoT d (x) = -eSf'(x) + o(e). 

where o is a continuous function with o(0) = 0. Observe that 5f > near rrii by the 
definition of m 8 . Plugging the above expansion in the ODE satisfied by T d , we see that: 

f + + f 
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We then take the limit as x \ rrii, so that e \ and — >• by the definition of mj. 
This leads to T' d (x) — > -1/2. 

Finally, we compute T^'(mj). By the ODE satisfied by and the smoothness of g, it 
follows that is different iable at any x > rrii. We then differentiate the ODE satisfied 
by T d , and use Taylor expansions. The result follows from direct calculation by sending 
x \ rrii. D 



3.5 Martingale version of the Brenier Theorem 

We next introduce a remarkable triple of dual variables corresponding to a smooth coupling 
function c. Recall the set D defined in (3.16) on which we have T d (x) = T u (x) = x, x G D, 
and the right-continuous inverse functions T^^T" 1 defined in Remark 3.8 (ii) and Remark 
3.9 (iii). 

The dynamic hedging component /i* is defined up to a constant, on each continuity interval, 
by: 

K = Cx( " T ; ) ~^ ( " Td) on D c , K = Ko T- 1 + c y (., .) - c y (T~\ .) on D. (3.18) 
The payoff function ip* is defined up to a constant on each continuity interval by: 

4>l = c y (T-\.)-KoT- 1 on D c , ^ = c y {T~\ .) - K o T' 1 on D. (3.19) 
The corresponding function is given by: 

V?*(x) = E F *[c{X,Y)-rb*{Y)\X = x] (3.20) 
= g(x)(c(x, .) - ^*) °r„(a;) + (1 - q(x))(c(x, .) - tp*) oT d (x), iGR. 

Finally, we define and ^* from (3.18)-(3.19) by imposing that 

the function c(.,T u ) — ip*{T u ) — [c(.,T d ) — yb^(T d )] — (T u — T d )h is continuous. (3.21) 

Theorem 3.13. Let fi, v be as in Assumptions 3.3 and 3.6. Assume further that(p+ e h l (n), 
ipt £ L 1 (z/) ; and that the partial derivative of the coupling function c xyy exists and c xyy > 
on R x R. Then: 

(i) (tp*,i>*, K) e v 2 , 

(ii) the strong duality holds for the martingale transportation problem, P* is a solution of 
P2{fJ>,v), and (</?*, tp*, /i*) zs a solution of D 2 (fJ,, v): 

J c(x,T*(x,dy))n(dx) = E ¥ *[c(X,Y)] = P 2 (/i,z/) = D 2 (n,v) = //(</?*) + 

Remark 3.14 (Mirror coupling: the right-monotone martingale transport plan). 

(i) Suppose that c xyy < 0. Then, the upper bound P 2 (/i,z/) is attained by the right-monotone 

martingale transport map 

F*(dx,dy) := Jx(dx)T^(x , dy) , 
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where T* is defined as in (3.15) with the pair of probability measures (fJ>,v): 

F p (x) := l-F^(-x), and F D {y) := 1 - F v {-y). 

To see this, we rewrite the optimal transportation problem equivalently with modified inputs: 

c(x,y):=c(-x,-y), #((-oo, x]) := fi([-x, oo)), &((-oo, y\) := v([-y, oo)), 

so that c xyy > 0, as required in Theorem 3.13. Note that the martingale constraint is 
preserved by the map (x,y) — > (—x, —y). 

(ii) Suppose that c xyy > 0. Then, the lower bound problem is explicitly solved by the right- 
monotone martingale transport plan. Indeed, it follows from the first part (i) of the present 
remark that: 

inf E F [c(X,Y)]=- sup E F [-c(X,Y)] =E K [c(X,Y)] = fc(x,T*(x,dy))n(dx). 

Remark 3.15. The martingale counterpart of the Spence-Mirrlees condition is c xyy > 0. 
We now argue that this condition is the natural requirement in the present setting. Indeed, 
the optimization problem is not affected by the modification of the coupling function from c 
to c(x, y) := c(x, y) + a(x) + b(y) + h(x)(y — x) for any a e b G L 1 (z/) 7 and h e L°. 

Since c xyy = c xyy , it follows that the condition c xyy > is stable for the above transformation 
of the coupling function. 

3.6 Comparison with Beiglbock and Juillet [4] 

The notion of left-monotone martingale transport was introduced by Beiglbock and Juillet 
[4], with an existence and uniqueness result, see Theorem 1.7 and Theorem 6.2. 

1. We first show that their conditions on the coupling function fall in the context of our 
Theorem 3.13: 

• The first class of couplings considered in [4] is of the form c(x, y) = h(y — x) for 
some differentiable function h whose derivative is strictly concave. Notice that 
this form of coupling essentially falls under our condition c xyy > 0. 

• The second class of couplings considered in [4] is of the form c(x,y) = ip(x)(f)(y) 
where ip is a non-negative decreasing function and a non-negative strict concave 
function. This class also essentially falls under our condition that c xyy > 0. 

2. The proof of [4] does not use the dual formulation of the martingale optimal transport 
problem. They rather extend the concept of cyclical monotonicity to the martingale 
context. As a consequence, [4] only provides an existence result and does not contain 
any explicit characterization of the maps (T^, T u ) and the optimal semi-static hedging 
strategy (</?*, ip*, h*). 
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Figure 1: Maps and T u built from two log-normal densities with variances 0.04 and 0.32. 
mi = 0.731. 

3. Our left-monotone martingale transport map T* coincides with the left-curtain coupling 
whose existence (and uniqueness) is stated in Theorem 4.18 of [4]. 

4. Our construction agrees with the example of two Log-normal distributions fio = 
e M{-ol/2,al) anc j Uo _ e Ar(-<j2/2,a2)^ ff 2 < ^ iU us trated in Figure 2 of [4]. By using 

our construction, we reproduce the left-monotone transference map in Figure 1. In- 
deed, in this case, Xq = — oo, SF has a unique local maximizer mi, which is then the 
global maximizer of SF, and X\ = oo. The left-monotone transport plan is explicitly 
obtained from our construction after Step 1, i.e. no further steps are needed in this 
case. 

Example 3.16. We provide an example where SF has two local maxima and the construction 
needs two steps. Let fi and v be defined by 

^1=7^(1,0.5) and v x (x) = - 2) +/\f(0.6, 0.1) + N(1A, 0.3)] . 

3 

Clearly \i and v have mean 1, and \i^v. We also immediately check that SF has two local 
maxima m\ = —0.15 and m<i = 0.72. Figure 2 below reports the maps T u and as obtained 
from our construction. 

3.7 Comparison with Hobson and Neuberger [32] 

Our Theorem 3.13 does not apply to the coupling function c(x,y) = \x — y\ considered 
by Hobson and Neuberger [32]. More importantly, the corresponding maps T^ N and T^ N 
introduced in [32] are both nondecreasing with T™(x) < x < T^ N (x) for all x G M. So 



17 



0.15 



4 - 



0.10 





-0.05 



-0.10 




-2 



-4 



Figure 2: 5F has two local maxima (left), and Td,T u corresponding to /ii, v\ (right). 

our solution (Td, T u ) is of a different nature and in contrast with the above (Tj* N , Tj* N ), our 
left-monotone martingale transport map T* does not depend on the nature of the coupling 
function c as long as c xyy > 0. 

However, by following the line of argument of the proof of Theorem 3.13, we may recover 
the solution of Hobson and Neuberger [32]. As a matter of fact, our method of proof is 
similar to that of [32], as the dual problem D2 is exactly the Lagrangian obtained by the 
penalization of the objective function by Lagrange multipliers. 

3.8 Some examples 

Example 3.17 (Variance swap). The coupling in this case is c(x,y) = In 2 (-) where ji and 
v have support in (0, 00). In particular, it satisfies the requirement of Theorem 3.13 that 
Cxyy > 0. Then, the optimal upper bound is given by 



where q is set to an arbitrary value on D. In Figure 3, we have plotted (p*,tjj* and h* with 
marginal distributions /i = e A ^" <T i //2,CT i) and u = e -^(-o-|/2,cr|) ^ ^2 _ < a 2 _ 32 y\? e 
recall that the corresponding maps T d ,T u are plotted in Figure 1. The expression for ip* is 




(3.22) 
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Figure 3: Superreplication strategy for a 2-period variance swap given two log- normal den- 
sities with variances 0.04 and 0.32. 



In particular, ip'l{x) 



for all x < mi . 



Example 3.18 (c(x,y) = — (f) P , p > 1, and \x,v have support in (0, oo)). This payoff 
function also satisfies the condition of Theorem 3.13 that c xyy > 0. The upper bound is 



X 



X 



fj,(dx). 



4 The n— Marginals Martingale Transport 

In this section, we provide a direct extension of our results to the martingale transportation 
problem under finitely many marginals constraint. Fix an integer n > 2, and let X = 
(X\, . . . ,X n ) be a vector of n random variables denoting the prices of some financial asset 
at dates t\ < . . . < t n . Consider the probability measures /i = (/ii, . . . ,/i n ) G (Vr) 71 with 
A*i zi ■ ■ ■ /4i in the convex order and 

J |£|/ij(cf£) < oo and J ^Hi(d^) = X , for all i = l,...,n. 

Similar to the two-marginals case, we introduce the set 

M n (fM) := {P G Vnifx) : X is a P- martingale}, 

where V n (fi) was defined in (2.6). In the present martingale version, we introduce the 
one-step ahead martingale transport maps defined by means of the n pairs of maps (T^, T*): 



Ti{x u 



W{*d + 1 D-{qi(x i )5 TUxi) + (1 - 9i)(xi)^(xi)) : 



(4.1) 
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where q t (0 := (f - Tj(0)/(^ - Tj)(0 for £ G A c , and (A,73,^)i=i,...,n-i are defined as 
in Subsection 3.3 with the pair (/ij,/i i+ i). 
The n— marginals martingale transport problem is defined by: 

P n M= sup E p [c(X)], 
PeA4„(M) 

where the map c : R n — >■ R is of the form 

for some upper semicontinuous functions c ! : 1 x 1 — > R with linear growth (or the 
condition (2.2)), i = \,...,n — \. 
The dual problem is defined by 

n 

D n (fi) : = inf V^(^), 

(u,h)(LV n ' 

where u = (tii, . . . , w n ) with components w J : R — )• R, and /i = (/ii, . . . , h n -i) with compo- 
nents hi : R* — >■ R, taken from the set of dual variables: 

V n :={(u,h): ( Ui )+ G L 1 ^), h t G L°(R*), and ©? =1 Ui + YJlZl hf > c}. 

Here, ©™ =1 w;(:r) = X^n^O^) and /if (x) = hi(xi, . . . ,Xi)(x i+ i - Xi). 

Similar to the two-marginals problems, the weak duality inequality P n (fi) < D n (fj) is 
obvious, and we shall obtain equality in the following result under convenient conditions. 

To derive the structure of the optimal hedging strategy, we shall consider the two-marginals 
(/ii, fii+i) problems with coupling functions c\ By Theorem 3.13, we have for i — 1, . . . , n — 1: 

PUm, := sup E p [c i {X,Y)}= inf .{(J>i(<p) + = Viiv*) + ft+^D) 

where i s defined as in (3.3) with c* substituted to c, and (ip*,ijj*, h*) G £>2 are defined as 
in (3.18)-(3.19)-(3.20) with d substituted to c and (T*, Tj) substituted to (T u , T d ). Finally, 
we define: 

u-(xi) := l {i < n }^*(xi) + l{i>i}^*_i(a;i), i = l,...,n, 
and m* := (uj, ...,<), /i* := (hi, . 

Theorem 4.1. Suppose fJ>i ^ . . . ^ fx n in convex order, with finite first moment, . . . , /i n _i 
/lave no atoms, and let Assumption 3.6 hold true for 5F = F^ i+1 — F^, for all 1 < i < n. 
Assume further that 

• c l have linear growth, that the cross derivatives c l xyy exist and satisfy c xyy > 0, 

• satisfy the integrability conditions (y?*) + G L 1 (/ij) ; G L 1 ( / Uj + i). 

Then, the strong duality holds, the transference map ¥* n (dx) = Hi(dx\) YYi=i ^*( a; i> ^j+i) ^ s 
optimal for the martingale transportation problem P n (fj), and (u*, h*) is optimal for the dual 
problem D n (fj), i.e. 

F* n eM n (n), (u*,h*)eV n , and E p "[c(X)]=P„(/i)=A t (/i) = Er=i^K*)- 
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Proof. Clearly, we have P* G .M„(/i), which provides the inequality E p «[c(X)] < P n {n). 
We next observe that (u*,h*) G T> n from our construction. Then D n (/j.) < Y2i<nf J ' i ( u i) = 
E p ™[c(X)]. The required result follows from the weak duality inequality P n (jj) < D n (/i). 



Remark 4.2. The optimal lower bound for a coupling function as in Theorem 4-1 is attained 
by the mirror solution introduced in Remark 3.14- 

Example 4.3 (Discrete monitoring variance swaps). This is a continuation of our Ex- 
ample 3.17. Suppose that fj, ± ^ . . . ■< /i n have support in (0, oo) with mean X , and let 
c(x u ...,x n ):=Zl 1 {^^tf- Then: 

Pn(j*) = J + [%(0(ln^) 2 + (l-^)(0(ln^) 2 ]^0- 

This optimal bound depends on all the marginals. The optimal lower bound is attained by 
our mirror solution, see Remark 4-2. 

Remark 4.4. In particular, their argument holds whenever c(x,y) which satisfies c(x,x) = 
= Cy(x,x), (x — y)c xy + c x > and our generalized Spence-Mirrlees condition c xyy > 0. 
Note that apart from the last condition, these requirements on c are not preserved by the 
transformation in Remark 3.15. 

Remark 4.5. In a related robust hedging problem, Hobson and Klimmek [31], derived an 
optimal upper bound for a derivative c(x±, . . . ,x n ) = Yl^i c °i x i, x i+i)- The difference with 
our problem above is that they are only given the marginal distribution jj, n for X n . See 
also Kahale [36]. We would like to emphasize that [31] assume the variance Kernel Co to 
satisfy the conditions c°(x,x) = c y (x,x) = 0, (x — y)c xy + c x > 0, together with our Spence- 
Mirrlees condition c xyy > 0. In the context of our problem with finitely many given marginals 
fj,i, . . . ,/i n , notice that, apart from the Spence-Mirrlees condition, none of these requirements 
are preserved by the transformation of Remark 3.15. 



5 Proof of the main results 

5.1 Construction of the left-monotone map 

This section is devoted to the proof of Theorem 3.10. We first motivate the definition of the 
maps T d and T u through the functions g and G. In this heuristic discussion, we consider the 
simple case of one single maximizer mi with 5F strictly increasing before mi, and we ignore 
the possible jumps of F v . 

The first observation about our construction is that for a point y G R, there are two 
alternatives: 
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• either y G (-00, mj; then P*[F G dy] = dF^y) + E[(l - g)(X)l {Td(x)G ^ } ] , and the 
requirement that Y v together with the decrease of Td imply that 

d(6FoT d ) = -(l- q )dF li ; 

in particular, in order for Td to be well-defined, it has to be valued in the domain of increase 
of 5F, 

• or y G (mi, 00), then P*[Y G dy] = E[q(X)l{ Tu ( X )ed y }], and the requirement that 
Y yU 2 together with the increase of T M imply that 

d(F u oT u ) = qdF^. 

Direct manipulation of these two equations implies that d5F(T d ) = —dF^ + dF v (T u ). Since 
Td{m\) = T u {mi) = mi, this implies that: 

F u (T u (x)) = F^x) + 5F(T d (x)), 

i.e. T u = g(.,T d ) as in (3.11), (3.12), and (3.14). 

Also, as a consequence of this relation, we see that the requirement T u (x) > x implies that 
SF(x) < 8F(T d (x)). Consequently, the choice of the break point mi as the maximizer of SF 
is necessary. 

We next substitute q and T u in the martingale condition: 

xdF^ = TvqdFp + TiidFp-qdF,,) = g(.,T d )d[F fl + 5F(T d )] — T d d5F(T d ). 
This implies that, for x > m\. 

F'^F^x) + y)dyj - T d d5F(T d ). 

Integrating from mi to x, and using the condition T d (mi) = mi, this provides: 

G(T d (x),x)-G(m 1 ,m 1 ) = 

where 

[K 1 ° ^(0 " iW^Z) + / F-\F,(x) + y)dy - / £d5F(Z) 

-00 J0 J(-oo,t] 

£dF„(0 - [ X £dF„(fl + / [g(x,0 - QdSF®, 

(-oo,F~ 1 oF fl (x)] J-00 J(~oo,t] 

in agreement with our definition of G A in (3.7) for A = (—00, mi]. We finally notice by direct 
computation that G(mi,mi) = 0, so that T d must satisfy the equation G(T d (x),x) = for 
all x >m\. 

Proof of Lemma 3.5 (i) Since SF is strictly increasing on A, we see that F v is strictly 
increasing in A. Therefore, for t < m < x, t £ A we have g(x,t) — t > g(t,t) — t = 0, 
implying that 1 1 — > G A (t,x) is strictly increasing in t on the set A. 
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We next verify that G A (m,x) > as long as SF(m) > SF(x). Denoting by d x the 
differential with respect to the x— variable, we compute by using the conditions on the set 
A that 



d x G A (m,x) = (F; 1 o F,{x) - x + [ dsgfaOUiOdSF^dF^x) 

V J(-oo,m] ' 

= (g(x,m)-x+ / dg{x,^jdF ll {x) > (g(x, m) - x)dF ll {x), 



since g(x,£) = (f>(5F(£)) where, for fixed x, the function y i-» <fi(y) := F~ 1 (5F(x) + y) is 
nondecreasing. Since F^ strictly increases at the right of m, and G A (m,m) = 0, this shows 
that G A (m, x) > as long as g(x,m) — x > or, equivalently, 8F{m) > SF(x). 

Then, in order to establish the existence and uniqueness of t A (x,m), it remains to verify 
that 



7 ( x ) : = G A (-oo,x) = / ZdF v (£)- £dF M (£)<0 for SF(x) < SF(m). 

J(-oo,F- 1 oF lx (x)] J-oo 

Let Xq := inf{x : SF(x) > 0}. Clearly, x < m < x, and 7 = on (—00,^0), j(r u ) = 0. 
Moreover, 7 is flat on Supp(F M ) c , where Supp(F M ) is a support of F^, and we see by direct 
differentiation that 7 is absolutely continuous with respect to \i with: 

dj(x) = (F- 1 oF t ,(x)-x)dF^(x), 

implying that dj < at the right of xo, by the (strict) convex-order property (/i ■< v) implied 
by the strict increase of SF on A. Furthermore, let x* be any possible local maximizer of 7. 
By the fact that 7 is flat off Supp(-F M ), we may assume that x* is either an interior point of 
Supp(F M ) or x* is a left accumulation point of Supp(F M ). In both cases, it follows from the 
first order condition that 

F- 1 (F,(x*)-)<x*<F- 1 (F,(x*)). 

If F~ x is continuous at the point F^(x*), then 5F(x*) = 0, and it follows from the definition 
of 7 that 



7 (x*) = / ZdSFfc) = - / (x*-Z)d8F(Z). 

J (— oo,x*] J (— 00, X*] 

By the (strict) convex-order property, this implies that j(x*) < 0. 

In the alternative case that F^ 1 jumps at the point F^x*), notice that F u is flat at the 
right of F^ 1 o F M (x*), and therefore the conclusion 7(2;*) < holds true in this case as 
well. Consequently, 7 < on (rr ,r M ). Since x > m > x , this provides the required strict 
inequality. 

(ii) Suppose x(m) < 00. Since 5F is strictly increasing on A, the inequalities (3.8) follow 
from the definition of x(m). 
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It remains to prove that 5F strictly increases in a right neighborhood of x(m) when- 
ever x{m) < oo. By definition, we have t A (x, m) > (SF)' 1 o SF(x) on (m,x(m)), and 
t A (x(m),m) < (5-F) -1 o SF(x(m)), where 5F" 1 denotes the inverse function of IacLSF. 
We denote h(x, m) := G ((8F)~ 1 o5F(x), x), and we compute that d x h(x, m) = [x — (5F)^ 1 o 
5F(x)]d5F(x). Since x > (5F)~ 1 o5F(x) whenever x > m, we see that h(.,m) decreases down 
from zero on the right neighborhood of x = m (confirming that t A (x,m) > (<5-F) _1 o SF(x) 
near m), and has the same maximum and minimum points as the function 5F. Since h must 
be increasing at a right neighborhood of x(m), it follows that 5F has the same property. 
D 

Proof of Theorem 3.10 (i) By construction, the probability measure P* satisfies the left- 
monotonicity property of Definition 3.2. In the rest of this proof, we verify that P* G 
-M 2 (/i, v). In particular, by the uniqueness result of Beiglbock and Juillet [4] (Theorem 1.5 
and Corollary 1.6), this would imply that P* is the unique left monotone transport plan. 

First, by the definition of P* in (3.15), X /z, and E P *[Y"|X] = X. It remains to 
verify that Y ~p„ v. We argue as in the beginning of Section 5.1 considering separately the 
following alternatives for any point t/Gl: 

Case 1: y = yd G D fl B corresponds to some point x such that yd = Td(x), and we see from 
the definition of P* that: 

P*[V G dy d ] = dF^Tdix)) - (1 - q)dF^x) and dF v (T u {x)) = qdF^. 

Then, F*[Y G dy d ] = dF^{y d ) - dF^x) + dF u (T u (x)). Since T u (x) = g(x,T d {x)), this 
provides P*[F G dy d ] = F v {dy) by direct substitution. 

Case 2: y = y u G D c corresponds to some x such that y u = T u (x), and we see from the 
definition of P* that: 

P*[Y G dy u ] = qdF^x) = d5F(T d (x)) + dF^. 
Using again the expression of T u in terms of T d , it follows that 

P*[Y G dy u \ = d(F u o T u (x) - F^x)) + dF^x) = dF u {x). 

Case 3: At a point of discontinuity of T u or T d , the above cases 1 and 2 are immediately adapted 
to account for the point mass. 

Case 4: In the remaining alternative y G D \ Bo, we observe that the function SF is flat near 
y, and there is no x ^ y such that T d {x) = y or T u {x) = y. Then, it follows from the 
definition of P* that: 

V.\Yedy]=dF li (y) = dF v (y). 
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(ii) Differentiating the integral equation defined by G A at a continuity point of T d , we see 
that: 

= -[F- 1 o F^x) - xjdF^x) + [g(x,T d (x)) - F' 1 o F^dF^x) 
+ [g(x,T d (x)) -T d (x)]dSF(T d (x)) 
= [g(x,T d (x)) - x\dF^x) + [g(x,T d (x)) - T d (x)]d5F(T d (x)). 

Since T u = g(.,T d ) this is the required ODE. The ODE for T u is obtained by using the 
relation T u = g (., T d ). O 

5.2 The optimal semi-static hedging strategy 

We start by following the same line of argument as in the proof of Theorem 2.2 in order 
to identify the semi-static hedging strategy introduced in (3.18-3.19-3.20). Our objective is 
then to construct a pair 

(</?*, e £>2 such that fifa) + z/O*) = E v *[c(X,Y)]. (5.1) 

This will provide equality in (3.5) with the optimality of P* for the optimal transportation 
problem P 2 and the optimality of (ip^,^, h*) for the dual problem D 2 . 
By the definition of the dual set V 2 , we observe that the requirement (5.1) is equivalent to 

V?*(X) + if>,(Y) + K{X)(Y - X) - c(X, Y) = 0, P* - a.s. for some function K,{5.2) 

and that the function is determined by: 

<p*(x) = m&xH(x, y), where H(x, y) := c(x, y) — ip*(y) — h*(x)(y — x), x,y<EM.. (5.3) 

The perfect replication property (5.2), is equivalent to: 

(p*(x) = q(x)(c(x,.)-ip„)oT u (x) + (l-q(x))(c(x,.)-^)oT d (x), (5.4) 
= (c(x,.)-^)oT(x)-(c(x )-^)oT d (x) 
1 ' (T u -T d )(x) 1 ; 

where we observe that we may choose arbitrarily on D. 

It remains to determine ip* by using the static superhedging condition (5.3). Since T u and 
T d are maximizers in (5.3), it follows from the first-order condition that 

ip'* ° T u (x) = c y (x, T u (x)) - K(x), V* ° T d (x) = c y (x, T d (x)) - K(x), x G D c , (5.6) 
and ip'^x) = c y (x,x) — h*(x) for x G D. (5.7) 

We now determine /i*. Differentiating (5.5), and using (5.6), we see that for x G D c : 

h , ^ = d f c(x, T u ) - c(x, T d ) | | T' u -T' d ip*(T u ) - ip*(T d ) 
dx I T u — T d J T u — T d T u — T d 
T>[c y (x,T d )-h]-T^[c y (x,T u )-h] 



T u — T d 
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Then, direct calculation leads to the expression of h'^ on D c reported in (3.18). Since T d and 
T u take values in D and D c , respectively, and h* is determined by the last two equations, we 
see that equation (5.6) determines on R. We finally observe that by (5.6) and (5.7), we 
have for x G D that ip'^(x) = c y (T^ 1 (x) , x) — h* oT^fi) = c y (x, x) — h*(x), which completes 
the definition of h*, up to an irrelevant constant, on D. 

5.3 Proof of Theorem 3.13 

Following the line of argument of the proof of Theorem 2.2, we see from the weak duality 
(3.5) that 

E p *[c(X,Y)] <P 2 (/^) < D 2 (ijl,u). 

Then, the proof of Theorem 3.13 is completed by the following result. 

Lemma 5.1. Let ji, v be as in Assumptions 3.3 and 3.6, and suppose that the payoff function 
c satisfies c xyy > 0. Then tp* © ip* + hf > c. 

Proof (i) We first verify that T u and satisfy the second order condition for a local 
maximum on D c . Differentiating (5.6), and using the expression of h'^ in (3.19), it follows 
from the condition c xyy > that, in the distribution sense, 

H yy {,T u )T' u = [c yy (.,T u )-^oT u ]Tl = Cx ^ T l ] ' C ^ Td) - c xy (.,T u ) < 

J-u — J-d 

H yy {.,T d )T< = [c yy (.,T d )-^oT d ]T d = ^ ^ ~ ^ Td) - c xy (., T d ) > 0, 

J-u — J-d 

on D c . By the nondecrease of T u and the nonincrease of T d , this implies that H yy (., T u ) < 
and H yy (.,T d ) < 0. 

(ii) We next show that y i — > H(.,y) is increasing before T d , and decreasing after T u . In 
particular, this implies that: 

(fJx) = max H(x,y) for all i6l. 

ye[T d (x),T u (x)] 

Set y := T u (x), let m« be the local maximum from which (T d ,T u )(x) is constructed, and 
consider an arbitrary y' = T u (x') > y for some x' > x. We only report the proof for the 
case x' G (rrij,Xj] for some j > i, the remaining cases are treated similarly. Recalling that 
H y (x,T u (x)) = 0, we decompose 

j 

H y (x,y') = H y (x,y') - H y (x,mj) +^2(A k + B k ), 

where 

A k := H y (x,m k ) - H y (x, z fe _i), B k := H y (x,x k _i) - i7 y (x,m fc _i A T u (x)). 
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We next compute from the expression of h* in (3.18) that: 

H y (x, y f ) - H y (x, mj ) = f [c yy (x,a - r(0]de 

Jmj 

< f [cyy(x,0-Cy y (T-\0,a}d? 

where the second inequality follows from the second order condition verified in (i). Similarly, 
we compute that 



A k = / [c yy (x,e)-r(e)]dc 

J Xfc_1 



/ ' Csyyfo <0, 

'Xk-1 Jx 

where we used again the second order condition verified in (i). Finally 



Bi 



= f [c yy (x,a-m')]de 

•/m fc _lVT u (rr) 

< f [c yy (x,a-c yy (T-\a,e)}de 

ry rTu^y') 

>' m fe-iVT u (x) J x 



A similar argument also shows that H y (x,y') < for y' < T d (x). 

(iii) We next show that H(.,T d ) = H(.,T U ). Denote 5H := H(.,T U ) — H(., T d ), and compute: 

6H' := c x (.,T u ) - c x (.,T d ) - (T u - T d )ti* 

+ [c y (.,T u ) - #(T U ) - K]T' U - [cy{.,T d ) - ^(T d ) - K]T' d 

in the distribution sense. By definition of ip* and ft,*, it follows that 5H' = at any 
continuity point. Since 5H is continuous by our construction, see (3.21), this shows that 
5H(x) = 5H(m,i) = 0, where m 8 is the local maximizer from which (T d ,T u )(x) is defined. 

(iv) We finally show that T u and T d are global maximizers of y i — > H(., y). Let x G D c , and 
denote by m the local maximizer from which T d (x) and T u (x) are constructed. For fixed 
T = T u {t) E (m,T u (x)) , it follows from similar calculations as in the previous step that 

d x {H(.,T u ) - H(.,T)} = c x (.,T u )-c x (.,T)-(T-T d )K 

trp _ T) f c x (.,T u )-c x (.,T) _ c x (.,T u ) -c x (.,T d ) \ 
{ u >\ T u -T ~ " ~ T u -T d ) 
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by the condition c xyy > 0. Then H(.,T U ) - H(.,T) = J t d x {H(.,T u ) - H(.,T)} > 0. 

By a similar calculation, we also show that H (x, Ta(x)^ — H(x,T) > for all T e 
(T d (x), m), thus completing the proof that T d and T u are global maximizers of y i — > H(., y). 

O 



6 Complement: Some examples of martingale mea- 
sures given marginals 

6.1 Local volatility model 

A first example of a martingale measure fitted to two marginal distributions /i tl and [i t2 , 
corresponding to the maturities ti < t 2 , is given by the Dupire local volatility model (in 
short LV) [23]. We first define an interpolation (^t)te[h,t 2 ] which does not violate the 
no-arbitrage condition, i.e. which obeys to the convex ordering condition. This can be 
achieved by introducing the implied Black-Scholes accumulated variances w(ti,K), defined 
by BS(K,w(ti,K)) = C(t u K) := f(^- K) + fi h (d£), where BS denotes the Black-Scholes 
formula 



with N the c.d.f. of the standard Normal distribution, and for t G [ti,t 2 ]: 

w{t, K) = l -^-uj{t x , K) + ^^w(t 2 , K), C(t, K) :=BS(X, zu(t, K)) = [ 

The Dupire LV model corresponding to this interpolation is defined by the SDE: 

o d t C(t,K) 



dX t = X t a loc (t,X t )dW t with a loc (t,KY 



d 2 K C(t,KY 



whenever <r loc is well-defined and induces a well-defined weak solution for the above SDE. 
In general, a\ oc is a measure with poor regularity. A rigorous adaptation of this solution, 
by convenient regularization of <7i oc , is provided by Hirsch and Roynette [28], resulting in a 
new proof of the Kellerer theorem [37]. 

A natural extension of such a LV model is given by the so-called local stochastic volatility 
model in which X t satisfies a non-linear McKean SDE ([26]): 

dX t = X t ^&2La t dW t 

with a t a (possibly multi-dimensional) Ito diffusion. Existence and uniqueness for such a 
non-linear SDE is not at all obvious and still open. 
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6.2 Local variance Gamma model 

A second example, introduced by P. Carr [9], which does not require the construction of 
a continuous-time implied volatility surface is given by the local variance Gamma model 
in which the process X t is defined as a time- homogeneous one-dimensional Ito diffusion X t 
subordinated by an independent Gamma process T t [9]: 

X t = X Ft 
dX t = a(X t )dW t , X = X 

The distribution of the Gamma process at time t is a Gamma distribution with density: 

cat* * i 

F{T t e ds} = wjTS sTW e ~ as > s > 

for some parameters t* — t 2 — h, a = l/t*. The Fokker-Planck PDE reads 

—<t[K) d K L{t 2 ,K) = t 2 -h 

from which we can deduce the local volatility function er(-) from call options valued uniquely 
at ti and t 2 . The Dupire infinitesimal calendar spread gets replaced by a discrete calendar 
spread. Similar to the previous example, a rigorous adaptation of this idea requires a 
regularization of the above function er(-) as in [28]. 

6.3 Local Levy's model 

As a last example, we review the local Levy model introduced by Carr, Geman Madan, and 
Yor [10]. The process X t is a compensated jump martingale 

dX t = I X t - (e x — 1) (m(dx, dt) — v{dx, dt)) , u(dx, dt) = a(t, X t )k(x)dxdt, 

where m(dx, dt) is the counting measure with compensator v. The analogue of the Dupire 
formula is 

d 2 y C(t,y)ya(t,y)ij(\n-)dy, 
with the double tail ip of the Levy measure k(x) given by 

/z px poo poo 

e x / k{u)dudx + l{ z >o} / e x k(u)dudx. 
-oo J— oo Jz Jx 
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